Newtonian N-body simulations are compatible with cosmological perturbation theory 
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Contrary to recent claims 1 1 1, Newtonian N-body simulations of collisionless Dark Matter in a ACDM back- 
ground are compatible with general relativity and are not invalidated by general relativistic effects at the linear 
level. This verdict is based on four facts. (1) The system of linearized Einstein equations and conservation laws 
is well-posed in the gauge invariant formulation and physically meaningful. (2) Comparing general relativity 
with its Newtonian approximation at a given order in perturbation theory is only meaningful at the level of 
observables. (3) The dynamics of observables describing a dust fluid in general relativity and its Newtonian 
approximation agree at the linear level. Any disagreement for observables on the lightcone are well-known, of 
which the most dominant is gravitational lensing. (4) Large fluctuations in the Hubble parameter contribute sig- 
nificantly only to gravitational lensing effects. Therefore, these fluctuations are not in conflict with Newtonian 
N-body simulations beyond what has already been carefully taken into account using ray tracing technology. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

There seems to be little doubt that the gravitational forma- 
tion of structures on scales deep inside the Hubble volume 
is accurately described by Newtonian gravity, with general- 
relativistic effects entering at subleading level in the pertur- 
bative description. Nevertheless it was claimed IH recently 
that effects in cosmological perturbation theory become dom- 
inant over 2nd order effects in the Newtonian approximation 
on scales larger than lOMpc. 

In greater detail the argument was based on the following. 
The evolution of cosmological perturbations was calculated in 
a particular coordinate system, called the Newtonian matter 
gauge (NM), in which the linear density contrast £nm and the 
peculiar velocity vnm coincide with the corresponding quan- 
tities £n = £nm, vn = vnm in the Newtonian (N) approxima- 
tion. It was found that fluctuations in the local Hubble param- 
eter S^ M = -1 - K NM /(3H), where K NM denotes the extrinsic 
curvature, are of order 0(Snm)- The significance of this obser- 
vation was evaluated by comparing c^ M to the size of second 
order corrections 8$ in the Newtonian approximation. For 
comoving scales k and redshifts z with 

^ M >4 2) and 4f«£ N <l, (1) 

it was argued that relativistic effects linear in cosmological 
fluctuations dominate over nonlinear Newtonian effects well 
within the domain of validity of perturbation theory. Based on 
this criterion, the authors of 1 1 ] found that linear cosmological 
perturbations on an Einstein-de Sitter background dominate 
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over Newtonian nonlinearities for scales k~ l > lOMpc during 
the redshift interval z e [0.4, 750], from which they concluded 
that Newtonian N-body simulations cannot be trusted on these 
scales during the specified redshift interval. 

Since by choice of gauge, 6^ = #nm> v n = vnm, while 
the fluctuations in the Hubble parameter are absent in the 
Newtonian approximation, and thus in N-body simulations, 
it is interesting to ask how these fluctuations become mani- 
fest in observables? It might be expected [ 1 ] that fluctuations 
in the Hubble parameter cause additional redshift space dis- 
tortions, because the observed redshift depends on K^m inte- 
grated along the line of sight between observer and source, in 
addition to the peculiar velocities of observer and source, v° M 
and v^ M , respectively. We find indeed that these fluctuations 
contribute considerably to redshift space distortions, however, 
only via gravitational lensing, which is well known as lensing 
magnification, e.g [2-8]. These lensing induced distortions 
are taken into account in N-body simulations using ray trac- 
ing technology l9lfT0l. 

Cosmological perturbation theory can be compared to its 
Newtonian approximation in a meaningful way only by com- 
paring observables. Observables are by definition gauge in- 
variant combinations of the perturbations. The set of observ- 
ables should be specified before a choice of gauge is imple- 
mented, after the gauge redundancies have been removed it is 
impossible to identify observables. However, once the observ- 
ables have been identified, gauge freedom is not sacrosanct 
and can either be removed or, equivalently, used to rewrite 
the theory using gauge invariant variables. Both procedures 
are perfectly valid. A comparison of observables using gauge 
invariant perturbations has been performed in |[TT1l and was 
criticized in (TJ as follows: "The initial conditions must be 
specified on a spatial Cauchy hypersurface, which in the con- 
text of cosmological perturbation theory corresponds to a par- 
ticular foliation of space-time, i.e., to a hypothetical observer 
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D = S h - 3<Hv L 
V = v L 



V = v NM 



FIG. 1 . D and V agree on different hypersurfaces (full Newton matter 
gauge and dashed longitudinal gauge) evaluated at the same coordi- 
nate values in their respective coordinate systems. Only the physical 
interpretation of D and V changes and on both hypersurfaces an hy- 
pothetical observer could determine D and V. 



who is able to determine physical quantities on a spatial hy- 
persurface. The relativistic-Newtonian correspondence mixes 
the quantities defined on different spatial hypersurfaces and 
thus no hypothetical observer in the Einsteinian world could 
actually determine these combined quantities." 

This is a misconception that requires immediate clarifica- 
tion. Although gauge invariant variables might have a sim- 
ple physical interpretation only in one specific gauge, they 
present observables in all other coordinate systems and, in 
particular, different hypersurfaces, as well. For instance, let 
D be a gauge invariant variable which reduces to the linear 
density contrast Ss in synchronous (S) gauge, measured by an 
observer at rest with respect to synchronous coordinates. Let 
V denote a gauge invariant variable which reduces to the pe- 
culiar velocity vl as measured by an observer at rest relative 
to the longitudinal coordinate system. This is the situation 
referred to "defined on different spatial hypersurfaces" in the 
above quote. However, (D,V) are defined in all possible co- 
ordinate systems and observers adopted to different and arbi- 
trary coordinate systems can measure D and V and will find 
the same numerical values for them. By construction, it is not 
the definition of gauge invariant variables that is tied to cer- 
tain hypersurfaces, but it is their physical interpretation. This 
is why any smart observer will construct the set of observables 
before choosing a particular gauge to measure them, because 
even if not all observables have a convenient physical interpre- 
tation in the observer's coordinate frame, they still resemble 
the only meaningful quantities for any other observer. The 
only reasonable academic debate between observers adopted 
to different coordinate systems is about the physical interpre- 
tation of the gauge invariant variables. 

As an example, consider an observer who is adopted to a 
longitudinal (L) coordinates system. This observer will in- 
terpret D as D L = £ L - 3*Hvl. Although the observer has 
a physical interpretation for the gauge dependent quantities 
vl), (s)he understands the necessity to construct gauge in- 
variant combinations involving (£ L , Vl), rather than assuming 
any other observer adopted to an arbitrary coordinate system 
would agree on the values of the gauge dependent quantities. 
The observer adopted to the longitudinal frame can measure 
v l) on a hypersurface, set up the initial value problem for 
an appropriate evolution equation involving (D,V) and finally 
solve for them, instead of (£l, vl), see Fig.[T] 

The plan for the rest of this paper is as follows. In Section|II| 
it is shown explicitly that the hypothetical observables (£n, vn) 
of linearized Newtonian gravity and (D,V) of linearized Ein- 



steinian gravity obey the same evolution equations on a flat 
ACDM background (in [ 1 ] only A = has been considered). 
This was expected since (D,V) reduce to (£nm,vnm) for ob- 
servers at rest in the Newtonian matter coordinate system. 

However, a real observer cannot measure (D,V) on a hy- 
persurface, because (s)he can only observe the light cone as- 
sociated with (her) him via light rays, which get affected by 
curvature perturbations. This induces a feedback of general 
relativistic effects on (D,V), although (D,V) still coincide with 
(£nm> v nm) on each hypersurface. Note that the status of (D,V) 
as observables is not challenged by the practical obstacles that 
prevent any observer from measuring them on the entire hy- 
persurface. 

Section [III] is devoted to investigate the impact of Hubble 
parameter fluctuations on the linear density contrast observed 
along a light cone. This is a well-known example highlight- 
ing the fact that S^ M contributes significantly only to gravita- 
tional lensing. The strongest additional redshift space distor- 
tions due to gravitational lensing GHU originate indeed from 
c^ M (as well as from intrinsic curvature, see Section III). The 
gauge invariant lensing term, of course, does not depend on 
the gauge the observer preferred. An observer adopted to the 
longitudinal gauge finds a negligible contribution from 6^ to 
gravitational lensing. 

Our main result is that Newtonian N-body simulations are 
in congruence with cosmological perturbation theory and are 
not threatened by relativistic effects at the linear level although 
relativistic effects can become significant. We give a pratical 
dictionary to use N-body simulation data to evaluate these cor- 
rections. 



II. CONGRUENCE OF LINEAR OBSERVABLES ON A 
HYPERSURFACE 

In this section we show that a pressureless fluid in a uni- 
verse with ACDM background geometry can be characterized 
by observables (D,V) in lst-order cosmological perturbation 
theory that obeys evolution equations identical to those gov- 
erning the evolution of (£n, vn). 

Restricting attention solely to scalar perturbations, the con- 
formal evolution equations for (0 N , £ N , vn) in the Newtonian 
approximation are given by 



A0 N = l<H 2 Q m 6 N 
6' N + Av N = , 

V' N + <Hv N = -0 N , 



(2a) 
(2b) 
(2c) 



with *H denoting the conformal expansion rate of the back- 
ground determined by Friedmann equations 



3<W 2 = 87rG(£> m + g A )a 4 , 
4nGa 2 g m = §£ m <W 2 = a 2 (<H 2 - <H') , 



(3a) 
(3b) 



where g m oc a~ 3 and Q m denotes the background matter den- 
sity relative to the critical density. 

The Newtonian perturbation variables (£n, vn) are defined 
by the dark matter density g m = g m (l + #n) and the peculiar 
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velocity Vn = Vvn- The triple (0n,£n,vn) constitutes the set 
of observables relevant for the discussion. 

The corresponding description in general relativity requires 
a background metric around which the geometry fluctuates. 
Restricting attention again solely to scalar metric fluctuations, 
(0, w, 0, h), the total metric field is give by 

ds 2 = a 2 [ - (1 + 20) dr dr + 2wj dr ® dx { + (4) 

+ [(1 - 2iff)6ij + 2hijW <g> dx j ] . 

The total pressureless source is encoded in T = g m U ® U, 
with g m = g m (l + 6) and U = (1 - 0, Vv)/a. Altogether the 
dynamical degrees of freedom (0, 0, w, h, S, v) are the gauge 
dependent metric, density and velocity perturbations. 

Observables can be constructed by the following gauge in- 
variant linear combinations: 

= + [(w - h')a]'/a , (5a) 

¥ = - <H(w - /*') , (5b) 

D = 8 - 3<H(v + w) , (5c) 

V = v + h! . (5d) 

The perturbed Einstein and conservation equations then yield 
evolution equations ifTTl for the gauge invariant quantities 
(0,¥,D,V): 

AO = I^OnD , (6a) 
D' + AV = , (6b) 

V + <WV = -O , (6c) 

where the background equations ([3]) and the linearized (0 /)- 
and (/j)-Einstein equations have been used: 

<H$ + = -l<H 2 Q m V , (7a) 
O = ¥ . (7b) 

Comparing ([6]) with the Newtonian approximation ^ it is ev- 
ident that the evolution equations are identical in form. In 
additions, (0, D, V) constitute the triple of observables rele- 
vant for this discussion. Of course, any linear combination 
of these gauge invariant variables constitutes an equally le- 
gitimate observable. The triple (0, D, V) is favored only to 
establish directly the correspondence between relativistic ob- 
servables and observables in the Newtonian approximation at 
the linear level. 

Note that (0, *F) has the same quasi-static dynamics as 0n, 
which allow us to qualify relativistic corrections to Newtonian 
observables, e.g. ^ below, as large or small in comparison to 
2nd order corrections in the Newtonian approximation ([T]). 

Let us emphasize again that there is no gauge Q in which 
simultaneously 

O = 0£, ¥ = 0£, D = %, V = Vg. 

Different observers simply assign different physical meaning 
to these gauge invariant variables. For instance, in longitudi- 
nal, synchronous and Newtonian matter gauge the following 
physical interpretations hold: 
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III. LINEAR OBSERVABLES ON THE LIGHTCONE 

A physical observer is in practice not able to measure (D,V) 
everywhere on any hypersurface. Instead, a physical observer 
can only learn about (D,V) by employing light rays travel- 
ing along (her) his respective light cone. As a consequence 
of such an observation campaign, (D,V) become subject to 
relativistic effects that have no Newtonian counterpart. The 
light rays will be gravitationally lensed and these lensing ef- 
fects will be attributed to (D,V). Since gravitational lenses 
are absent in the Newtonian approximation, the dictionary 
(O, D, V) <-> (0n, vn) is challenged. 

As an example, consider the linear density fluctuation 
A(#i, z) at the observed redshift z and direction -n on the sky 
or, equivalently, in the direction n of the incoming light ray 
at the observer's space-time position (to, x<j), within the rela- 
tivistic framework of 1st order cosmological perturbation the- 
ory. It is given by the gauge invariant expression Q 

A(n, z) = D - j-d 2 Y - — [ dA-^^A Q {<$> + Y) + 
r s Jo r 

+ (I£ + ^)( t -* v+ .^*H (8) 

i 2 r s 

<H r s Jo 

where all functions are evaluated along the unperturbed 
light cone x = x - nr(z), r = r - r(z), with r(z) = 
f Q z dz f /(H(z')ao), the unperturbed affine parameter X - r and 
S denotes the source. A^ is the angular part of the Laplacian 
in spherical coordinates. Detailed derivations of ^ can be for 
instance found in HE1E1. 

The Newtonian approximation of gravity is void of the light 
cone concept. Moreover, the gravitational potential couples 
only to massive bodies, in particular, there is no coupling to 
light rays. In other words, light rays cannot probe 0n- The 
assumption that there is a light cone embedded in the back- 
ground cosmology, albeit an artificial point of view, allows 
to obtain the first and second term of ^ in the Newtonian 
approximation, A^(n,z) - £n - *H -1 ^vn. The second term 
is known as Kaiser-effect lfT2l and is the dominant redshift 
space distortion for small redshifts. Moreover, let us assume 
that light couples to 0n such that its bending around the Sun 
conforms to actual observations. Including the lensing contri- 
bution, 

A N (n, z) = S N - <H~ 2 d 2 r v N - — [ dA r -^A Q( p N . (9) 

rs Jo r 

Depending on scale, redshift and redshift binning, the lens- 
ing contribution can be the leading redshift space distortion, 
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which can even dominate over c) N Q[8) for sufficiently distant 
sources. 

For a more transparent treatment, let us define a Newtonian 
observable An(w, z) through the following replacements in the 
relativistic quantity ((SI: 



tures as follows: 



A N (w,z) = A(n,z)\ 



lD-xJ N ,V->v N ,O=Y->0 N 



(10) 



Using the results from the last section it follows that 
An(w, z) = A(n, z). As a consequence, N-body simulations can 
be used directly to extract relativistic observables when scalar 
dust fluctuations on a ACDM background are considered at 
the linear level. 

Let us comment on why fluctuations c^ M in the Hubble pa- 
rameter K = -3H, 



rH_ SK/3 

6 =-^Ta 



-1 



(if/ +<H(p+ 1/3 A(w -/*')) 



(11) 



are, in fact, strongly contributing only to the lensing term and, 
therefore, were identified correctly in [lj as a major source 
of redshift space distortions. Clearly, large fluctuations in the 
Hubble parameter do not imply that Newtonian N-body sim- 
ulations cannot be trusted. Instead it implies that either ^ 
or ( fT0| ) (or, better, an expression including nonlinear effects) 
should be used to compare numerical experiments based on 
the Newtonian approximation to observations, which was well 
known GHE) before the work Q. 

We could have argued based on gauge invariance alone 
(3j Q that A(n, z) is constructed from extrinsic curvature 8 H , 
intrinsic curvature = 4Aif//a 2 , anisotropic extrinsic curva- 
ture Aij = a(didj/A - 8ij/3)A(w - h'), and the divergence of 
the observer's coordinate acceleration a = Vln[a(l + 0)] in 
such a way that arguing about the size of 8 H adapted to vari- 
ous gauges is meaningless. In certain gauges 8 H might qualify 
as large (NM and S), while in others it qualifies as small (L), 
but this does not matter at all. 

Any change in 8 H induced by a transition between coordi- 
nate systems is compensated for by corresponding changes in 
8 R = Aifr, 8 a = A0 and 8 A = A(w - h'). This can be checked 
explicitly for the gauge invariant lensing term A(0+ x F), which 
can be deconstructed into the various gauge dependent curva- 



A(0 + Y) = 8 R + 8 a + & 



A' 



(12) 



Note that 8 H contains the anisotropic extrinsic curvature 8 A . 
In the Newtonian matter gauge, 8^ M becomes large just be- 
cause c^ M is large, while in the longitudinal gauge 8 A = and 
0l, <Al are quasi-static, meaning they basically remain at their 
initially small values (except close to neutron stars and black 
holes, see [Q21 [141 for more details). As a consequence, 8^ 
is negligible. Since 0l = ® and = all other terms in 
^ involving *F and O, as well as their conformal time deriva- 
tives remain much smaller than the first three terms in ([8]). It 
can be shown that whenever c^ M contributes to these less rel- 
evant terms, its c^ M component is either compensated for by 
another c^ M contribution or rendered harmless by A -1 . 

The discussion outlined for A(n, z) applies quite generally 
to any observable. A constructive algorithm for an arbitrary 
observable is the following: (/) construct the general rela- 
tivistic gauge invariant observable and express it in terms 
of (D,V) and O. (ii) Use the quasi- static evolution of O to 
determine which contributions qualify as large. If one can 
identify large contributions that are not reflected in the corre- 
sponding Newtonian observable, then these contributions will 
have a genuine relativistic origin. (Hi) Employ the dictionary 
(D, V, O) — > (£n, vn, 0n) to extract relativistic observables us- 
ing Newtonian N-body simulations. This has been worked out 
in much greater in detail in (T5l[T6ll . 



IV. CONCLUSION 

In summary, we have shown that fluctuations in the Hub- 
ble parameter do not give rise to new, so far overlooked red- 
shift space distortions. They contribute to redshift space dis- 
tortions, however only to the well-known lensing magnifica- 
tion. Therefore contrary to the claims of Q, Newtonian N- 
body simulations are the appropriate numerical experiments 
to extract information on relativistic observables at least in 1st 
order perturbation theory. 
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